Examining coupled-channel effects in radiative charmonium transitions 
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Coupled-channel effects due to coupling of charmonia to the charmed and anticharmed mesons are of cur- 
rent interest in heavy quarkonium physics. However, the effects have not been unambiguously established. 
In this Letter, a clean method is proposed in order to examine the coupled-channel effects in charmonium 
transitions. We show that the hindered Ml radiative transitions from the 2P to IP charmonia are suitable 
for this purpose. We suggest to measure one or more of the ratios V{h'^ -> XcJ7)/r(Xcj ~^ Xcj""") and 
r(Xc.7 ~^ hc'y)/T{x'cj — ^ XcJi""), for which highly nontrivial and parameter-free predictions are given. The 
picture can also be tested using both unquenched and quenched lattice calculations. 
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Thanks to various experiments world-wide, our knowledge 
of the physics of heavy quarkonium has been greatly enriched 
in the last decade. New charmonium(-like) states, includ- 
ing the he and x'c2 as well as the so-called XYZ states, 
were observed. Most of the XYZ states are above the open- 
charm thresholds, and do not fit the expectations from the 
quark model. Hence, it is of current interest and high im- 
portance to investigate the coupled-channel effects, originat- 
ing from the coupling of cc to charmed-meson-anticharmed- 
meson channels, in charmonium physics. So far, these ef- 
fects have not been established unambiguously, though ev- 
idences exist, see, e.g^ Refs. Ill-Ill for transitions between 
charmonia and Refs. iQ. Iltl for spectroscopy. This Letter is 
devoted to a clean way of examining the coupled-channel 
effects. For this purpose, we propose to measure one or 
more of the ratios T{h[. Xcjj)/^{x'cj ~^ Xc/tt*^) and 
^ix'cj ~^ ^c7)/r(Xcj — ^ Xcj7r°), for which highly nontriv- 
ial predictions will be made. 

At the hadronic level, one may consider the coupling 
of heavy quarkonium states to open-flavor meson and anti- 
mesons using effective Lagrangians, and take into account 
the coupled-channel effects by calculating intermediate heavy 
meson loops. Because the difference between the mass of a 
heavy quarkonium and heavy meson-antimeson thresholds is 
small, the intermediate heavy mesons are nonrelativistic with 
velocity v 1. Based on this observation, a nonrelativis- 
tic effective field theory (NREFT) was proposed |[lllll3]. It 
was found that the transitions between two P-wave charmonia 
with the emission of one pion are completely dominated by 
the coupled-channel effects with an enhancement of 1 /v^ 
in the decay amplitude in contrast to the multipole contribu- 
tion HI]. Here we find that the hindered Ml transitions from 
a 2P to a IP charmonium is also dominated by the coupled- 
channel effects. Because these two types of transitions can- 
not be directly connected to each other, only when both of 
them are dominated by the coupled-channel effects, nontrivial 
predictions can be made in the framework of NREFT. Con- 
sequently, the measurements of hindered Ml transitions from 
P-wave charmonia provides a good opportunity for examin- 



ing the coupled-channel effects in charmonium transitions. 

There are several nice features of the hindered Ml tran- 
sitions of P-wave charmonia for investigating the coupled- 
channel effects of open-charm mesons: 
1) First of all, these transitions are expected to be dominated 
by the coupled-channel effects. On one hand, in quark models, 
the decay amplitude for an Ml transition between two heavy 
quarkonia is proportional to the overlap of the wave functions 
of the initial and final quarkonia (see, e.g. 181]), 
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with ipi(f) being the wave function of the initial (final) heavy 
quarkonium, and the photon energy in the rest frame of the 
initial particle. For transitions between a 2P and a IP state, if 
the charmonia are purely cc states, the overlap is nonzero only 
because of small relativistic corrections. The statement can be 
made model-independently using the potential nonrelativistic 
QCD ^ — the leading contribution vanishes Holl . Hence, the 
transition amplitude would start from E^Vc/rric, where Vc and 
rric are the charm quark velocity and mass, respectively, and 
the factor 1 /nic accounts for the spin-flip. Indeed, the tran- 
sition rates are very small in quark model calculations — at 
the largest of the order 1 keV On the other hand, because 
the leading coupling of a P-wave charmonium to a charmed- 
meson and anticharmed-meson pair is in an S'-wave, the de- 
cay amplitude through intermediate charmed-meson loops as 
shown in Fig.[T](a) scales as 
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where and [v^Y account for the nonrelativistic loop 
measure and three nonrelativistic propagators, respectively, 
and E^ comes from the P-wave coupling of the photon to 
charmed mesons. The factor of l/?7ic is again due to the 
spin-flip (more details will be given below Eq. (O). In addi- 
tion, the amplitude is proportional to the electric charge e and 
the product of coupling constants of the IP and 2P states to 
the charmed mesons. Yet, there is no suppression analogous 
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FIG. 1: Possible triangle (a) and two-point (b) loops for the radia- 
tive transitions, (c) and (d) are two typical two-loop diagrams. The 
double, solid, wavy and dashed lines represent charmonia, charmed 
mesons, photons, and pion, respectively. 



to because the initial and final charmonia, though 

having different principal quantum numbers, do not couple to 
each other directly. Instead, they couple through intermedi- 
ate charmed mesons, and there is no similar suppression for 
such couplings. Note that v in Eq. (O will not approach 
even when the charmonium mass overlaps with the charmed- 
meson threshold, since it should be understood as the average 
of two velocities corresponding to the two cuts in the three- 
point loop. 

2) The triangle hadronic loops involved in the transitions are 
convergent in the nonrelativistic framework. Therefore, we do 
not need to introduce a counterterm. On the contrary, similar 
loops for the Ml transitions of 5'-wave charmonia are diver- 
gent so that similar statements cannot be made there. This 
is another nice example of the important role of such dia- 
grams in hadron physics, see e.g. the classical work on neu- 
tral pion photoproduction off nucleons lITltl or the more re- 
cent investigation of the large isospin violation in the decay 
77(1405/1475) ^ Stt ifll . 

3) Because the leading coupling of charmed and anticharmed 
mesons to the P-wave charmonium is in an S'-wave, there 
is no derivative in such vertices. Hence the two-point loop 
with the four-particle contact term XcjD^*^ D'-*'>j, as shown 
in Fig.[T](b), is not related to that in Fig.[T](a) by gauge sym- 
metry. It can be treated separately. Being gauge-invariant by 
itself, the four-particle contact term should contain the electro- 
magnetic field strength F^'^ for the photon. Hence, the corre- 
sponding vertex is proportional to the external momentum of 
the photon. The amplitude for the diagram scales as 
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One sees that it is two orders higher in the meson velocity 
counting than the diagram Fig. [T] (a), and hence can be ne- 
glected at leading one-loop order 

4) As will be shown later, the two-loop diagrams (c) and (d) 
are also suppressed compared with (a). 

The coupling of the P-wave charmonia to the charmed and 
anticharmed mesons is described by the Lagrangian 113 1 



■ 91 



Tr [x^'Haa'Ha] + h.c. 



(4) 



where gi is the coupling constant of the IP charmonium 
states (g[ will be used for the 2P states). Ha = Va ■ a + Pa 
and Ha — — Va • <? + are fields annihilating charmed and 
anticharmed mesons, respectively, with a the Pauli matrices 
and a the light flavor index. The two-component notation in- 
troduced in Ref. |fT^ is used here, which is convenient for 
processes with negligible recoil effect (less than 1% for the 
processes considered in this Letter). The P-wave charmonia 
are collected in the spin-multiplet 
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The magnetic coupling of the photon to heavy mesons is de- 
scribed by the Lagrangian lfl4lll5ll 
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Tr 
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(6) 

where B'^ = g^jkQij^j j-j^g magnetic field, Q = 
diag{2/3, —1/3, —1/3} is the light quark charge matrix, and 
Q' is the heavy quark electric charge (in units of e). The first 
term describes the nonperturbative physics of the light quarks, 
while the second term is for the magnetic coupling of the 
heavy mesons and hence is proportional to l/niQ. Although 
the photon can also couple to the heavy mesons through gaug- 
ing the kinetic energy term, this vertex does not contribute to 
the magnetic transitions. 



K 7XcO 

K 7Xc2 



[D,D',D'],[D',D,D],[D\D\D*] 
[D,D*,D], [D*, 0,0*1 [0*,D*,0] 
[D,D*,D*], [D*,7?*,D*] 
[D*,D*,D], [D*,D*,D*] 



TABLE I: Possible loops contributing to each transition. The charge- 
conjugated ones and the flavor labels are not shown for simplicity. 



Because the Xcj states are easier to be detected than the he, 
and x'c2 has been observed, we will calculate the decay widths 
of the hindered Ml transitions h'^ — >■ jXcj and x'c2 ~^ 7^c- 
Results for the other hindered Ml transitions of the P-wave 
charmonia can be easily obtained using the same method. De- 
noting the charmed meson connecting the initial charmonium 
and the photon as Ml, the one connecting two charmonia as 
M2, and the other as M3, we specify the triangle loops by 
[A/1, M2, A/3]. Considering both the pseudoscalar and vec- 
tor charmed mesons, possible loops for these transitions are 
listed in Table|T](see also Ref. jit] for details). 

The decay amplitude for each transition can be expressed 
in terms of the scalar three-point loop function 
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where P and q are the momenta of the initial particle and the photon, respectively, mi{i = 1,2, 3) are the masses of the particles 
Mi in the loop. The analytic expression can be found in Refs. |[ll[l6t]. The amplitude for the transition x'c2 ~^ l^c reads 
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where the loop function has been written as 
I{q, Ml, M2, M3). The charge-conjugated channels 
are taken into account.' The amplitudes for the other 
transitions can be obtained similarly. 

Since the spin direction of the c or c quark should be flipped 
in the Ml transitions, the decay amplitude should vanish in the 
heavy quark limit. It is nonzero only because of the 0{m'^^) 
spin symmetry breaking effect. One easily sees there must be 
nonvanishing contributions from the second term in the La- 
grangian Eq. (|6]l. In fact, the first term, to be called /3-term in 
the following, also contributes at the same order though m^^ 
is not explicit in the amplitude. Let us look at the decay am- 
plitude given in Eq. (O. The /3-term contribution would van- 
ish were spin symmetry a good symmetry, i.e., different loops 
proportional to /3 cancel each other exactly if the hyperfine 
splitting between vector and pseudoscalar charmed mesons 
Muj^j — Md* ^ is tuned to zero. The surviving part is due to 
the nonvanishing hyperfine splitting which is of order m^'. 

Because the expansion parameter in the NREFT v ~ 0.4 
is not small, the results should have sizeable uncertainties. 
This can be seen by analyzing the power counting of the de- 
cay amplitudes for certain two-loop diagrams. In Ref. flii, it 
is argued that vertex corrections due to pion-exchange is sup- 
pressed, so that the largest two-loop contribution comes from 
diagrams shown in Fig.[r|(c) and (d). The diagram (c) contains 
four nonrelativistic charmed meson propagators and one rela- 
tivistic pion propagator. Each momentum is of order Mjjv, so 
that each propagator scales as 1/w^ in the velocity counting. 
The photon vertex comes from gauging the charmed-meson- 
pion axial coupling, hence it contributes a factor of g / F^, with 
g and F.^ being the axial coupling constant and pion decay 
constant, respectively. The charmed-meson-pion axial cou- 
pling is in a P wave. Because this is the only P-wave vertex 
in the diagram, it should scales as the photon momentum, and 
the vertex is proportional to E^g/ F^^. Therefore, the decay 
amplitude for the diagram shown in Fig.[T](c) scales as 
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' They were not considered in Ref. flll^ ItI [TtIi . Hence all the loop am- 
plitudes therein should be doubled, and the decay widths from the loops 
should be multiplied by 4. All ratios remain innocent. We thank T. Mehen 
and D.-L. Yang for pointing out this. 



where the factor l/(47r)2 appears because there is one more 
loop than in the one-loop case, and the chiral symmetry break- 
ing scale is = 47rF7r- In order to compare with Eq. (|2|i, 
a factor of M]^ is introduced, with M^ being the charmed- 
meson mass, to make the whole scaling have the same dimen- 
sion as that in Eq. (|2]i. The value of the axial coupling constant 
g = 0.6 can be determined from r(D* Dtt) UJi, J^. Nu- 
merically, one has gMu/A^ ~ 1. The diagram (d) has the 
same scaling as (c). This can be seen easily because the one 
more propagator in (d) is balanced by two more P-wave ver- 
tices. Therefore, the two-loop diagrams are effectively sup- 
pressed compared with Fig.[T](a) by a factor of t; ~ 0.4. 

In numerical calculations, we use the central values of all 
measured masses lfl9ll . The value of /? is not precisely known. 
Here, we take the value = 276 MeV determined with 
iric = 1.5 GeV in Ref. il411 . In fact, the precise value of f3 is 
not important. A change of from 276 to 376 MeV only 
causes a change in the decay width of less than 10%. Hence, 
the decay width for the x'c2 ^ 7^c is 



r(xc2 ihc) 



= (10.7±4.3)-^^i^keV, 



(10) 



where a 40% uncertainty has been assigned to account for 
higher order effects. Because the XcJ and he are below the 
open charm thresholds, the coupling constant gi cannot be 
measured directly through the decays of the P-wave charmo- 
nia. Similarly, since the only established 2P charmonium x'c2 
is below the D*D* threshold, g[ is also not known yet. In 
order to obtain an order-of-magnitude estimate of the branch- 
ing fraction for the decay x'c2 ~^ l^c, we take a model value 
gi ~ —4 GeV^'^/^ 2 p^j. ^/ ^ resort to quark model 
calculations of the decay widths of x'cj ™d ^c- From the 
results in the nonrelativistic potential model in Ref. j^, we 
get g[ = 0.5. ..1.3 GeV^^^^. From the results in the Cornell 
coupled-channel model 0], we get g[ = 0.7. ..1.2 GeV^^/^ 
Hence, one may take g[ ~ 1 GsV~^^^ as an estimate. Thus, 
we get the estimate T{x'c2 ^ l^c) = C(170 kcV), which is 
much larger than the quark model prediction 1 .3 keV iK] . This 
is consistent with the analysis made above that the transitions 
are dominated by the coupled-channel effects. The width of 
the x'c2 was measured to be (24 ± 6) MeV lfl9ll . Hence, the 



■ The value of gi as defined in Eq. (4) is twice of that in flsll . 
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FIG. 2: (a) Decay widths of the Ml transitions h'^ — !■ "fXcj- (b) 
Parameter-free ratios r{h'^ 7Xcj)/r(Xc2 '"""xci)- The dot- 
dashed, sohd and dashed lines are for h'^ decays with XcO, Xci and 
Xc2 in the final state, respectively. 



branching fraction for the Ml transition x'c2 ~^ 7^c is 
B{x'c2 ^ ih,) = 0{1 X 10-2). 



(11) 



A future observation with comparable branching fraction 
would strongly indicate the dominance of coupled-channel ef- 
fects in the transition. Similar predictions can be made for the 
other Ml transitions of P-wave charmonia. The results for the 
widths of the transitions h'^ IXcJ are shown in Fig.|2](a) as 
a function of the unknown mass of the h'^. They should be un- 
derstood to have an uncertainty of about 40%. One sees that 
the partial widths of these transitions are of the same order as 
r(Xc2 ^ l^c)^ and are typically of 0(100) keV if using the 
same model estimates of gi and 5ias above, much larger than 
the results in the potential model 181]. 

More interestingly, nontrivial parameter-free predictions 
can be made for the ratios of the partial widths of these hin- 
dered Ml transitions to those of the transitions between two 
P-wave charmonia with emission of one pion. In Ref. [jlj], it 
is shown that the latter transitions are also dominated by the 
charmed meson loops. The one-loop and largest two-loo p di - 
agrams scale as q^A/v^ iH and {q^A/v'^){MDE^/Al) U^, 
respectively, where q^[ET^) is the pion momentum (energy), 
and the charged and neutral charmed meson difference A de- 
scribes the isospin breaking. One sees that the one-loop di- 
agrams dominate over the two-loop ones. The decay ampli- 
tudes and decay widths for the single-pion transitions have 
been calculated in Ref. III. Here we only compare the Ml 
transitions with the decay x'c2 ~^ Xci7r°. The width was pre- 
dicted in Ref. as (0.29±0.10)(gigi GcV)^ kcV, see foot- 
note [U The right panel of Fig. |2] shows the parameter-free 
predictions of the ratios T{h'^ 7Xc,/)/r(Xc2 ^ ^°Xci) 
as a function of the mass of the h'^. Such predictions can 
only be made when both processes are loop-dominated be- 
cause only in this case the decay amplitudes are proportional 
to the same product of coupling constants .91.9^. Were they 
multipole-dominated, an unknown matrix element of gluon 
operators would be involved in the transition x'cz ^ '"'"Xci so 
that the process cannot be directly related to the hindered Ml 
transitions. Even though the uncertainty of these predictions 



is sizeable, they are markedly different from potential model 
calculations. 

In this Letter, we argue that the hindered Ml transitions of 
P-wave charmonium are dominated by the coupled-channel 
effects. The conclusion is supported by numerical calcula- 
tions. With a reasonable estimate of the unknown coupling 
constants, the results turn out to be much larger than those 
obtained in the quark model. Parameter-free predictions are 
made for ratios of partial widths of two completely different 
types of charmonium transitions: the hindered Ml transitions 
and single-pion transitions of the P-wave charmonia. The P- 
wave charmonia considered here are assumed to be cc states 
so that they are organized as in Eq. (|5]l. If their coupling to the 
charmed mesons becomes resonant, renormalization is nec- 
essary (see, e.g. Refs. ll20l 21]). In that case, the resulting 



widths for the hindered Ml transitions should still be much 
larger than the results were the transitions not dominated by 
the coupled-channel effects, since a resonant coupling tends 
to enhance the widths further. Experimental efforts on mea- 
suring the transitions suggested here are needed towards un- 
derstanding the coupled-channel effects in the charmonium 
transitions. In fact, the existing technologies of lattice cal- 
culations of the heavy quarkonia radiative transitions 1122 - 
2411 are well ready to test the picture presented in this Letter: 
if the hindered Ml transitions of the P-wave charmonia are 
dominated by the coupled-channel effects, the results of sim- 
ulations with dynamical light quarks should be significantly 
larger than those in the quenched approximation. 

We further notice that the hindered Ml transitions of the 5- 
wave charmonia are not well suited for studying the coupled- 
channel effects. This is because for these transitions the 
charmed meson loops are divergent. While the divergence of 
the triangle diagrams scales as 0(w") in the velocity count- 
ing, the finite part scales as 0{v). The divergence must be 
absorbed by a counterterm. However, the counterterm cannot 
be determined from elsewhere. Hence, parameter-free predic- 
tions analogous to those made in this Letter is not possible. 
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